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Abstract — Instability of population 1 (X = 0.7, Z = 0.02) massive stars against radial 
oscillations during the post-main-sequence gravitational contraction of the helium core is in- 
vestigated. Initial stellar masses are in the range GSM© < Mzams < 90Mq. In hydrodynamic 
computations of self-exciting stellar oscillations we assumed that energy transfer in the enve- 
lope of the pulsating star is due to radiative heat conduction and convection. The convective 
heat transfer was treated in the framework of the theory of time-dependent turbulent convec- 
tion. During evolutionary expansion of outer layers after hydrogen exhaustion in the stellar 
core the star is shown to be unstable against radial oscillations while its effective temperature 
is Teff > 6700K for Mzams = SSM© and T^s > 7200K for Mzams = 90Mq. Pulsational insta- 
bility is due to the K-mechanism in helium ionization zones and at lower effective temperature 
oscillations decay because of significantly increasing convection. The upper limit of the period 
of radial pulsations on this stage of evolution does not exceed ~ 200 day. Radial oscillations of 
the hypergiant resume during evolutionary contraction of outer layers when the effective tem- 
perature is Teff > 7300K for Mzams = SSMq and T^s > 7600K for Mzams = 90Mq. Initially 
radial oscillations are due to instability of the first overtone and transition to fundamental 
mode pulsations takes place at higher effective temperatures {Tes > 7700K for Mzams = Q5Mq 
and Teff > 8200K for Mzams = QOMq). The upper limit of the period of radial oscillations of 
evolving blueward yellow hypergiants does not exceed ~ 130 day. Thus, yellow hypergiants are 
stable against radial stellar pulsations during the major part of their evolutionary stage. 
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1 INTRODUCTION 



More than a half-century ago Feast and Thackeray (1956) pointed out the significant scatter 
in luminosities (up to 7 magnitudes) of supergiants of the Large Magellanic Cloud. The most 
luminous stars of this sample of spectral types F and G the authors of this work named as super- 
supergiants. Later the name supcr-supcrgiant was changed to a hypergiant and the luminosity 
class of such objects was denoted as la"*" or laO (dc Jager 1980). The high luminosity and the low 
surface gravity favour the strong stellar wind, so that as one more criterion of the hypergiant 
one should consider the widely broadened emission and absorption lines. Evolutionary status 
of yellow hypergiants are thought to be that they are the post-main-sequence massive stars 
(i.e. with initial masses Mzams > 4OM0) undergoing gravitational contraction of the helium 
core (Meynet et al. 1994; Stothers and Chin 1996; Langer et al. 2007). 

Yellow hypergiants are variable objects but the common point of view on nature of their 
variability does not exist yet. Stothers and Chin (1996, 2001) showed that outer layers of 
hypergiants might be dynamically unstable because the adiabatic exponent becomes smaller 
its critical value (Fi < 4/3) due to partial ionization of helium and the large contribution of 
radiation pressure. Outer layers of the hypergiant can be unstable also due to convection in 
the hydrogen and helium ionization zones. Dc Jager and Nieuwenhuijzen(1997) and de Jager 
et al. (2001) estimated that the gradient of the turbulent pressure in subphotospheric layers 
might reduce the effective gravity almost to zero. However this conclusion was not confirmed 
by Stothers (2003) who showed that the turbulent pressure is insignificant and this effect can 
be neglected. 

Together with irregular enhancements of the stellar wind some yellow hypergiants show 
cyclic light variations. The most remarkable example is p Cas with light amplitude of a few 
tenths of a magnitude and the period of order of several hundred days (Arellano Ferro 1985; 
Sheffer and Lambert 1986; Zsoldos and Percy 1991; Lobel et al. 1994, 1998). Most of these 
authors mentioned the fact that cyclic light variations of p Cas are very similar to those of pul- 
sating stars. Moreover, near-IR CO line profiles demonstrate splitting similar to that observed 
in pulsating supergiants of intermediate spectral types such as RV Tau and R CrB variables 
(Gorlova et al. 2006). At the same time there seems to be a contradiction with pulsation 
hypothesis because light variations with period of ^ 500 day cannot be explained in terms of 
radial stellar pulsations for typical masses and radii of yellow hypergiants (Sheffer and Lambert 
1986). Here one should be noted that no studies of pulsational instability of yellow hyper- 
giants have been done yet, so that disussion of their photometric variability in terms of stellar 
pulsations remains hypothetical. 

As shown in our previous paper (Fadeyev 2011) the population I stars {Z > 0.02) with 
initial mass Mzams > 6OM0 become unstable against radial oscillations during the main- 



sequence evolutionary stage and they pulsate at least up to hydrogen exhaustion in the stellar 
core. The goal of the present study is to consider the pulsation properties of massive stars during 
the later evolutionary stage of gravitational contraction of the helium core and to determine 
the role of radial stellar pulsations in variability of yellow hypergiants. 



2 THE METHOD OF COMPUTATIONS 

The methods used to compute the stellar evolution and self-exciting radial stellar oscilla- 
tions were described in our previous papers (Fadeyev 2007, 2008a, 2008b, 2010, 2011), so that 
below we outline the most important improvements employed in hydrodynamic calculations of 
the present work. In studies of radial pulsations of Wolf-Rayet stars (Fadeyev 2007, 2008a, 
2008b) and LBV-stars (Fadeyev 2010) the hydrodynamic computations were done in assump- 
tion of radiation transfer without convection. Bearing in mind complete ionization of hydrogen 
and helium in outer layers of these stars such an approach seems to be appropriate. In cal- 
culations of nonlinear pulsations of massive main-sequense stars with effective temperatures 
Tes > 2 X lO^K (Fadeyev 2011) the energy transport was assumed to be due to radiation and 
convection. The region of convective heat transfer encompasses the zone of the iron Z-bump 
(T ~ 2 X lO^K) and convection was treated in the framework of the local steady-state mixing 
length theory (Bohm-Vitense 1958). However in outer layers (i.e. in hydrogen and helium 
ionization zones) of yellow hypergiants with effective temperatures Tcs < lO^K the significant 
fraction of energy is transported by convection and interaction of convective elements with 
pulsation motions should be taken into account because of the large pulsation amplitude. To 
this end in the present study we employ the equations of time-dependent turbulent convection 
obtained by Kuhfuss (1986) from Navier-Stokes equations written for spherical geometry. 

For spherically-symmetric motions the Lagrangean mass coordinate Mr and the radius r 
are related by continuity equation dM^ — Awr^pdr, where p is the gas density. Equation of 
motion describing the radial gas flow velocity U as a function of time t is 

^ - i^,2 g(^ + ^t + Pv) _ GM^ 

dt dMr r2 ' ^ ' 

where P is the total thermodynamic (gas plus radiation) pressure, Pt is the turbulent pressure, 
Py is the turbulent viscous pressure, G is the Newtonian graviational constant. The energy 
conservation equation is 

«(^ + ^.)+(p^p.^p,)^m=_iv.F. (2) 



dt dt \p J p 

where E is the specific internal energy (the sum of the translational, excitation, ionization 
and radiation energies per gram of material), Et is the root-mean specific turbulent energy, 
F = Fiad + Fc + Ft is the total flux, that is the sum of the radiative, convective (enthalpy) and 
turbulent (kinetic energy of convective elements) fluxes. 



Additionally we write the equation for conservation of the turbulent energy 



(3) 



dt 

where Dt is dissipation of turbulent energy due to molecular viscosity and is the source or 
sink of the turbulent energy due to boyancy forces. In the isotropic medium the turbulent 
pressure Pt, the mean turbulent energy Et and the root-mean speed of convective elements Vc 
are related by 

Pt = ^Et = pvl (4) 

Expressions for other quantities in equations (II]) - ([3]) can be found in papers by Wuchterl and 
Feuchtinger (1998), Olivier and Wood (2005), Smolec and Moskalik (2008). 

The numerical solution of equations (II]) - (IH]) was done with difference equations of the 
second-order accuracy in both the spatial coordinate and the time t. The equation of 
motion ([1]) was solved with the explicit finite-difference method, whereas equations ([2]) and 
([3]) were treated implicitly using the Crank-Nicholson scheme (see, e.g., Richtmyer and Morton 
1967). Thus, at each step of integration with respect to time t the temperature T and the specific 
turbulent energy E^ are determined for all Lagrangean mass zones of the hydrodynamic model 
from iterative solution of linearized difference energy equations. 

For the inner boundary condition we employed the assumption of the rigid sphere with the 
time-independent flux: 

where tq and Lq are the radius and luminosity at the inner boundary. The gas temperature 
at the inner boundary is T < lO^K and therefore the thermonuclear energy generation within 
the pulsating envelope can be neglected. In hydrodynamic calculations the luminosity at the 
inner boundary Lq was treated as the initial parameter of the hydrodynamic model determined 
from computations of stellar evolution. It should be noted that for the sake of stability of 
the numerical solution the inner boundary condition of the energy equation was written in 
assumption of radiation heat conduction. Therefore, the inner boundary of the hydro dynamical 
model is always below the iron Z-bump convection zone (T ~ 2 x lO^K). For models considered 
in the present study the radius of the inner boundary is ro ^ 3 x 10~^i?, where R is the 
equilibrium radius of the outer boundary. 

At the evolutionary stage of the yellow hypergiant the stellar mass is about a half of its initial 
mass and in the outer layers of the star there is the negative gradient of the mean molecular 
weight. The spatial inhomogeneity of the chemical composition is due to diminution of the both 
stellar mass and mass of convective core during the main-sequence hydrogen burning. However 
the mass of outer pulsating layers of the hypergiant is enough small in comparison with stellar 
mass (Menv < 10~^M), so that effects of inhomogeneity of chemical composition with respect 
to spatial coordinate remain always negligible. 



3 RESULTS OF COMPUTATIONS 



3.1 Stellar evolution 

Yellow hypergiants are on the short evolutionary stage of gravitational contraction of the 
helium core during of which the star evolves along the loop in the Hertzsprung-Russell (HR) 
diagram and its effective temperature becomes as low as 4 x lO^K < T^s ^ 5 x lO'^K. The parts 
of the evolutionary tracks of stars with effective temperatures T^s < lO^K and initial masses 
A^^zAMS = 65Mq and QOM© are shown in Fig. [I](a). The time of evolution along the shown 
parts of the tracks ranges from ^ 10^ yr (Mzams = 90Mq) to ^ 2 x 10"^ yr (Mzams = Q5Mq). 

During evolutionary expansion of outer layers when the stellar effective temperature is T^s = 
lO^K (the initial point of the tracks in Fig. [It^a)) the principal energy source is gravitational 
contraction of the helium core. For example, the ratio of central values of the thermonuclear 
energy generation rate Ec (these are mostly reactions of the 3-a process) and the gravitational 
energy release ranges from ac = ed {—TcdSc/dt) ^ 7 x 10"'^ for Mzams = 65M0 to ~ 0.05 
for Mzams = QOM©. Here Tc and Sc are the temperature and specific entropy in the center 
of the star. This ratio gradually increases with time and at the turning point of the track is 
Oc ~ 10^ independently of Mzams- At ending points of the tracks shown in Fig. [T](a) this ratio 
is etc ^ 300 for Mzams = SSMq and W for Mzams = QOM©. 

Gradual decrease of the role of gravitational contraction in the total stellar energy produc- 
tion leads to deceleration of the evolutionary movement along the track. This effect is illustrated 
in Figim^b), where for tracks in Fig[T](a) we give the plots of the rate of effective temperature 
change dhiTcs/dt. For example, at Teg ~ 7- lO'^K the rate of the redward evolution in the HR 
diagram about five times of that in the opposite direction at the same effective temperature. 
Thus, the probability to observe the yellow hypergiant with increasing effective temperature is 
significantly higher than that to observe the hypergiant with expanding outer layers. Neverthe- 
less in our study we consider the pulsation properties of hypergiants with effective temperatures 
Tes < lO^K independently of direction of their evolution in the HR diagram. 

3.2 Pulsations of yellow hypergiants 

The study of self-exciting stellar oscillations implies solution of the Cauchy problem for 
equations ([I]) - ([3]) with initial conditions corresponding to the hydrostatic and thermal (V-F = 
0) equilibrium. To this end we used some stellar models of preliminary computed evolutionary 
sequnces. Stability of the star against radial oscillations is determined from integration of the 
equations of hydrodynamics with respect to time t. In the case of pulsational instability the 
solution describes the exponetially growing oscillation amplitude with following transition to 
the limit cycle oscillations. The main quantities characterizing the limit cycle are the period H 



and the amplitude of the radial displacement of the outer boundary AR expressed in units of 
the equilibrium stellar radius R. 

In Fig. |2] the plots of AR/R and 11 as a function of T^s are shown for the evolutionary 
sequence of models with initial mass Mzams = 65Mq. As in the HR diagram the effective 
temperature increases from right to left. The dashed lines correspond to the evolutionary 
expansion of outer layers of the star, that is to decrease of the effective temperature, whereas 
solid lines correspond to their contraction. 

Radial oscillations of the hypergiant at the initial point of the evolutionary track (Teg = 
lO^K) are driven by the layers of the iron Z-bump (T ~ 2 x lO^K) and by helium ionization 
zones locating closer to the surface. The role of the iron Z-bump in pulsational instability 
decreases with decreasing effective temperature because of diminution of the relative radius 
r/R and the amplitude of radial displacement Ar in these layers. That is why the evolutionary 
decrease of the effective temperature is accompanied by decreasing pulsation amplitude AR/R. 
This conclusion is illustrated in Fig. |3] where for two hydro dynamical models with effective 
temperatures T^s = lO^K and T^s = 6800K we show the radial dependencies of the mechanical 
work over a closed thermodynamic cycle: 



Here AMj is the mass of the j-th Lagrangean interval, Pj is the sum of the thermodynamic 
pressure and turbulent pressure, Vj is the specific volume. 



is the normalizing coefficient, is the number of Lagrangean zones in the hydro dynamical 
model. Depending on the stellar structure the hydrodynamic computations were done with 
500 < < 1000. 

As seen in Fig. [3] the evolutionary decrease of effective temperature from lO^K to 6800K is 
accompanied by decrease of the relative radius of the layer with maximum mechanical work Wj 
at the iron Z-bump from r/R ^ 0.3 to r/R ^ 0.15. At the same time the contribution of the 
iron Z-bump zone into the positive mechanical work decreases from 0.43 to 0.1, so that near 
the boundary of the pulsational instability radial pulsations are driven mostly in the helium 
ionization zones. Here we have to note that the plots of the mechanical work Wj in Fig. |3] are 
shown as a function of the equilibrium relative radius r/R, whereas the mechanical work of a 
group of mass zones is determined from integration with respect to mass coordinate Mr. 

Evolutionary expansion of the hypergiant outer layers is accompanied by increase of convec- 
tion in the helium ionization zones, so that for Tcs ^ 6700K the radial oscillations decay due to 
large decrease of the radiative flux in ionization zones. Recommencement of radial oscillations 




(6) 




(7) 



becomes possible after the turn of the evolutionary track when the effective temperature rises 
up to Tcff > 7300K. Pulsational properties of the hypergiant during evolutionary contraction 
of outer layers are illustrated in Fig. HI where radial dependencies of the mechanical work over 
a closed thermodynamic cycle are shown for hydrodynamical models with T^s = 7400K and 
lO^K. 

The star of effective temperature T^s = 7400K with contracting outer layers locates near 
the boundary of the pulsational instability region. In comparison with the hypergiant of the 
same effective temperature but on the preceding evolutionary stage of expanding outer layers 
(the luminosity and the surface mass fraction of helium are L = 1.033 x IO^Lq and Yg = 0.74) 
the star has a somewhat higher luminosity and perceptibly higher surface abundance of helium: 
L = 1.085 X IO^Lq, Yg = 0.88. The higher helium abundance leads to a larger opacity of stellar 
material and to stronger convection in the helium ionization zones, so that the boundary of 
the pulsation instability moves to somewhat higher effective temperatures in comparison with 
hypergiants undergoing expansion of outer layers. 

Pulsational instability of the hypergiant during evolutionary contraction of its outer lay- 
ers is due to the K-mechanism in the helium ionization zones. The radial dependence of the 
mechanical work over a closed thermodynamic cycle for the model near the boundary of in- 
stability region is shown in Fig. |4] by the dashed line. Two prominent positive maxima on the 
plot correspond to the zones of the first ionization and second ionization of helim, whereas the 
contribution of the iron Z-bump zone into the positive mechanical work does not exceed 6%. 
During further evolution with increasing effective temperature the role of the iron Z-bump zone 
remains insignificant and at Tes = lO^K its contribution into the total positive mechanical work 
remain less than 8%. 

Of greatest interest in comparison with observations is the pulsation period 11. During 
the evolutionary expansion of outer layers in the hypergiant with initial mass Mzams = 65Mq 
the radial pulsations are due to instability of the fundamental mode and as seen in Fig. Mjo) 
the period gradually increases with decreasing Teg up to the instability boundary. During 
movement along the track in the opposite direction the radial oscillations resume at effective 
temperature T^s = 7400K. For T^s < 7700K the star pulsates in the first overtone and the period 
of the contracting hypergiant is significantly shorter than that of the expanding hypergiant with 
the same effective temperature. Transition of radial pulsations from the first overtone to the 
fundamental mode takes place at effective temperature Tes ~ 7700K and is revealed by a jump 
in plots of the surface amplitude AR/R and period 11 (see Fig. [2]). 

To compare the hypergiant pulsation properties at evolutionary stages of expansion and 
contraction we give in Fig. [5] the plots of the gas flow velocity at the outer boundary U and 
the light in V band my for two hydrodynamical models with initial mass Mzams = Q5Mq and 
effective temperature T^g = 7500K. The signiflcant difference in amplitudes of velocity and 



light curves is due to the fact that evolutionary contracting hypergiant pulsates in the first 
obertone, whereas during evolutionary expansion the hypergiant pulsates in the fundamental 
mode. However at effective temperatures Teg > 7700K the difference in pulsational properties 
becomes insignificant (see Fig. |6]) because on both evolutionary stages the hypergiant pulsates 
in the fundamental mode. 

Above we discussed results of hydrodynamic calculations for stars with initial mass Mzams = 
65Mq. Hypergiants with larger initial mass have a higher helium abundance in their outer 
layers because of the more massive convective core and the stronger stellar wind during the 
main-sequence evolutionary stage. Therefore in more massive hypergiants with Tgg < lO^K the 
convection in helium ionization zones is stronger and the boundary of the pulsation instability 
region corresponds to higher effective temperatures. For example, in stars with initial mass 
^ZAMS = 9OM0 the radial pulsations cease during evolutionary expansion of outer layers at 
effective temperature Tes ~ 7200K and resume during evolutionary contraction at Tcs ~ 7600K. 

In Fig. [7] the amplitude of the radial displacement of the outer boundary AR/R and the 
pulsation period 11 are shown as a function of effective temperature for hypergiants with initial 
mass Mzams = QOMq. With approaching the boundary of pulsational instability during the 
evolutionary expansion of outer layers the mode of radial oscillations changes and hypergiants 
with effective temperature T^s < 7800K pulsate in the first overtone. Therefore, the maximum 
period of radially pulsating hypergiant with initial mass Mzams = QOMq does not exceed 
150 day. During movement along the track in the opposite direction the radial oscillations 
resume in the first overtone at effective temperature Tes ~ 7700K and at Tes ^ 8400K the star 
oscillates in the fundamental mode. 

The main properties of some hydrodynamic models of radially pulsating hypergiants with 
initial masses Mzams = SSM© and 90Mq are summarized in the table. All models are repre- 
sented in four sequences depending on the initial stellar mass Mzams and direction of movement 
in the HR diagram. The stellar mass M, the luminosity L and the surface fractional mass abun- 
dance of hydrogen Xg approximately do not change within each sequence and therefore are given 
only for the first model of the sequence. The luminosity L and effective temperature T^s cor- 
respond to the initial hydrostatic equilibrium, 11 and Q are the pulsation period and pulsation 
constant in days, k is the order of the pulsation mode {k = for the fundamental mode and 
/c = 1 for the first overtone). The rate of period change n/II = dlnU/dt was evaluated from 
the second-order numerical differentiation with respect to age of the stellar evolution model 
used as initial conditions in hydrodynamic calculations. Thus, to evaluate the rate of period 
change we had to compute two additional hydrodynamical models. In the table we give only 
most reliable estimates of this quantitiy. The amplitude of the radial displacement of the outer 
boundary AR expressed in units of the initial equilibrium radius R and the amplitude of the 
radial velocity at the outer boundary AU are given in last two columns of the table. 



4 Conclusions 



Results of hydrodynamic calculations described in the present paper allow us to conclude 
that radial oscillations is rather a short-term phenomenon in the life of yellow hyper giants. For 
example, hypergiants with initial masses 65Mq < AfzAMS < 90Mq and effective temperatures 
^eflf < lO^K are stable against radial oscillations during ^ 95% of the time needed to move 
along the loop of the evolutionary track. The stability against radial oscillations is due to the 
strong convection in helium ionization zones which suppresses the K-mechanism of pulsation 
instability, whereas the zone of the iron Z-bump locates to deep in the stellar envelope and 
does not play any perceptible role. Therefore, radial pulsations of yellow hypergiants take place 
only in the beginning and in the end of this evolutionary stage. 

Estimates of the upper limit of period of radially pulsating hypergiants obtained in the 
present study can be applied to interpretaion of some variable stars. For example, semi-regular 
light variations of p Cas with period of 300 day < H < 500 day definitely cannot be explained 
in terms of radial stellar oscillations. Another example is the pulsating variable star V1427 Aql 
with uncertain evolutionary status. In some studies this variable is discussed as the post-AGB 
stars, whereas in other works as a yellow hypergiant (see, for discussion, Le Coroller et al. 
2003). Multicolor photometric observations of V1427 Aql carried out during eight years reveal 
the presence of light variations with periods of 130 and 200 day, the star moving blueward in 
the HR diagram (Archipova et al. 2009). Results of our computations allow us to conclude 
that V1427 Aql is rather the post-AGB star because assumption on yellow hypergiant leads 
to following contradictions. First, the mean effective temperature of this star is Teg — 6750K 
(Kipper 2008) and corresponds to the region of pulsational stability. Second, the period of hght 
changes ~ 200 day significantly exceeds its upper limit for radially pulsating yellow hypergiants 
evolving blueward in the HR diagram. 

Discussed above results of hydrodynamic calculations deal with massive stars with effective 
temperatures Tea < lO^K. In the present study we carried out hydrodynamic computations 
for some models with higher effective temperatures corresponding to the both earlier and later 
stages of stellar evolution. In the first case pulsational instability allows us to conclude that 
radial oscillations of main-sequence stars with initial mass Mzams > 60Mq (Fadeyev 2011) exist 
during the post-main-sequence stage and decay only at effective temperature Teg ~ 7000K due 
to convection. In the second case radial oscillations arise at the final stage of the yellow 
hypergiant and do not decay up to exhaustion of helium in the stellar core. With increasing 
effective temperature the hypergiant becomes the LBV-star radial pulsations of which are 
responsible for the microvariability (Fadeyev 2010) and then transforms into the Wolf-Rayet 
star which is also pulsationally unstable (Fadeyev 2007, 2008a, 2008b). 



REFERENCES 

1. A. Arellano Ferro, MNRAS 216, 571 (1985). 

2. V.P. Arkhipova, V.F. Esipov, N.P. Ikonnikova, G.V. Komissarova, A.M. Tatarnikov, B.F. 
Yudin, Pis'ma Astron. Zh., 35, 846 (2009) [Astron. Lett., 35, 764 (2009)]. 

3. E. Bohm-Vitense, Zeitshrift fiir Astrophys. 46, 108 (1958). 

4. H. Le CoroUer, A. Lebre, D. Gillet, et al., Astron. Astrophys. 400, 613 (2003). 

5. Yu.A. Fadeyev, Pis'ma Astron. Zh., 33, 775 (2007) [Astron. Lett. 33, 692 (2007)]. 

6. Yu.A. Fadeyev, Astron. Zh. 85, 716 (2008a) [Astron. Rep. 54, 645 (2008)]. 

7. Yu.A. Fadeyev, Pis'ma Astron. Zh., 34, 854 (2008b) [Astron. Lett. 34, 772 (2008)]. 

8. Yu.A. Fadeyev, Pis'ma Astron. Zh., 36, 380 (2010) [Astron. Lett. 36, 362 (2010)]. 

9. Yu.A. Fadeyev, Pis'ma Astron. Zh., 37, 13 (2011) [Astron. Lett. 37, 11 (2011)]. 

10. M. W. Feast and A. D. Thackeray, MNRAS 116, 41 (1956). 

11. N. Gorlova, A. Lobcl, A.J. Burgasser, et al.), Astropliys.J 651, 1130 (2006). 

12. C. de Jager, The Brightest Stars (Reidel, Dordrecht, 1980; Mir, Moscow, 1984). 

13. C. de Jager, Astron. Astrophys. Rev. 8, 145 (1998). 

14. C. de Jager and H. Nieuwenhuijzen, MNRAS 290, L50 (1997). 

15. C. de Jager, A. Lobel, H. Nieuwenhuijzen, et al., MNRAS 327, 452 (2001). 

16. T. Kipper, Baltic Astron. 17, 87 (2008). 

17. R. Kuhfuss, Astron. Astrophys. 160, 116 (1986). 

18. N. Langer, C.A. Norman, A. de Koter, et al., Astron. Astrophys. 475, L19 (2007). 

19. A. Lobel, C. de Jager, H. Nieuwenhuijzen, et al, Astron. Astrophys 291, 226 (1994). 

20. A. Lobel, G. Israehan, C. de Jager, et al., Astron. Astrophys. 330, 659 (1998). 

21. G. Meynet, A. Maeder, G. Schaller, et al., Astron. Astrophys. Suppl. 103, 97 (1994). 

22. E.A. Olivier and P.R. Wood, MNRAS 362, 1396 (2005). 

23. R.D. Richtmyer and K.W. Morton, Difference Methods for Initial-Value Problems (2nd 
ed., Interscience, 1967; Mir, Moscow, 1972). 



24. Y. Sheffer and D. L. Lambert, PASP 98, 914 (1986). 

25. R. Smolec and P. Moskalik, Acta Astron. 58, 193 (2008). 

26. R.B. Stothers, Astrophys. J. 589, 960 (2003). 

27. R.B. Stothers and C.-W. Chin, Astrophys. J. 468, 842 (1996). 

28. R.B. Stothers and C.-W. Chin, Astrophys. J. 560, 934 (2001). 

29. E. Zsoldos and J. B. Percy, Astron. Astrophys. 246, 441 (1991). 

30. G. Wuchterl and M.U. Feuchtinger, Astron. Astrophys. 340, 419 



Hydrodynamic models of yellow hypergiants 



M'7 A l\/rO / Mr^ 


M/Mn. 

■1 V-L 1 -J. V-L f . J 


LlLr^ 


X 


eii ) 


n 


O 


k 


fi/n 


AR/R 


AU 






10*' 






(\(\\ 










kill / s 


65 


32.9 


1.03 


0.24 


10 


63 


0.0587 





3.6 


0.65 


140 










9 


82 


0.0551 





3.6 


0.28 


79 










8 


115 


0.0541 





2.4 


0.33 


84 










7 


168 


0.0525 





1.3 


0.33 


77 










6.8 


180 


0.0516 







0.28 


60 




28.6 


1.08 


0.10 


7.5 


110 


0.0380 


1 


-0.44 


0.08 


25 










8 


124 


0.0523 





-0.60 


0.37 


92 










9 


87 


0.0524 





-0.60 


0.31 


86 










10 


64 


0.0525 







0.28 


87 


90 


42.2 


1.52 


0.16 


10 


74 


0.0584 







0.42 


no 










9 


100 


0.0559 





2.5 


0.32 


91 










8 


143 


0.0561 







0.39 


101 










7.5 


124 


0.0402 


1 


1.5 


0.09 


27 




38.4 


1.52 


0.08 


7.7 


123 


0.0413 


1 




0.07 


19 










8 


111 


0.0421 


1 


-0.44 


0.14 


49 










9 


103 


0.0553 





-0.63 


0.36 


101 










10 


76 


0.0563 







0.31 


100 



FIGURE CAPTIONS 



Fig. 1. (a) - Parts of the evolutionary tracks of stars with initial mass Mzams = QOM© (solid 
line) and Mzams = 65Mq (dashed line) during the post-main-sequence core contraction. 
Arrows indicate the direction of evolution along the track, (b) - The rate of effective 
temperature change dlnTes/dt along the track. 

Fig. 2. The amplitude of the radial displacement of the outer boundary AR/R (a) and the 
period of radial pulsations 11 (b) as a function of the effective temperature Tgg for models 
of the evolutionary sequence Mzams = SSMq Dashed lines correspond to the evolutionary 
expansion of outer layers of the star and solid lines correspond to their contraction. 

Fig. 3. The normalized mechanical work over a closed thermodynamic cycle Wj as a function 
of the relative equilibrium radius r/RaX the evolutionary stage of expansion of outer 
layers of the star Mzams = GSM© at effective temperature T^k = lO^K (solid line) and 
Teff = 6800K (dashed line). 

Fig. 4. Same as in Fig. [3] but for the evolutionary stage of contraction of outer layers of the star 
with Teff = 7400K (dashed line) and Teg = lO^K (sohd line). 

Fig. 5. The gas flow velocity at the outer boundary U (a) and the light in V band my (b) of 
hydrodynamic models of the hypergiant with initial mass Mzams = 65Mq and effective 
temperature Tgff = 7500K. Solid and dashed lines correspond to evolutionary stages of 
contraction and expansion, respectively. 

Fig. 6. Same as in Fig. [5] but for hydrodynamical models of the hypergiant with effective tem- 
perature Teff = 8000K. 

Fig. 7. Same as in Fig. [2] but for hypergiants with initial mass Mzams = OOM©. 




Figure 1: (a) - Parts of the evolutionary tracks of stars with initial mass Mzams = 90Mq (solid 
line) Mzams — Q5Mq (dashed line) during the post-main-sequence core contraction. Arrows 
indicate the direction of evolution along the track, (b) - The rate of effective temperature 
change d In Tes/dt along the track. 




Figure 2: The amplitude of the radial displacement of the outer boundary AR/R (a) and the 
period of radial pulsations 11 (b) as a function of the effective temperature Tes for models of the 
evolutionary sequence Mzams = 65Mq Dashed lines correspond to the evolutionary expansion 
of outer layers of the star and solid lines correspond to their contraction. 
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Figure 3: The normalized mechanical work over a closed thermodynamic cycle Wj as a function 
of the relative equilibrium radius r/R at the evolutionary stage of expansion of outer layers of 
the star Mzams = OSM© at effective temperature Tes — lO^K (sohd line) and Tes — 6800K 
(dashed line). 




Figure 4: Same as in Fig. [3] but for the evolutionary stage of contraction of outer layers of the 
star with T^s = 7400K (dashed line) and Tgg = lO^K (solid line). 




Figure 5: The gas flow velocity at the outer boundary U (a) and the hght in V band my (b) of 
hydrodynamic models of the hypergiant with initial mass Mzams = 65Mq and effective tem- 
perature Tes — 7500K. Solid and dashed hues correspond to evolutionary stages of contraction 
and expansion, respectively. 




Figure 6: Same as in Fig. |5]but for hydrodynamical models of the hypergiant with effective 
temperature Tes = 8000K. 



